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Abstract 

Repetitive tasks in industrial works may contribute to health prob¬ 
lems among operators, such as musculo-skeletal disorders, in part 
due to insufficient control of muscle fatigue. In this paper, a pre¬ 
dictive model of fatigue is proposed for repetitive push/pull opera¬ 
tions. Assumptions generally accepted in the literature are first ex- 
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plicitly set in this framework. Then, an earlier static fatigue model 
is recalled and extended to quasi-static situations. Specifically, the 
maximal torque that can be generated at a joint is not considered as 
constant, but instead varies over time accordingly to the operator’s 
changing posture. The fatigue model is implemented with this new 
consideration and evaluated in a simulation of push/pull operation. 
Keywords '.Predictive model; Fatigue; Push/pull operations; Repeti¬ 
tive tasks; simulation 


1 Introduction 


A better understanding on how the human body deals with muscle fatigue 
would allow better estimates of working abilities and enhance production. 
Such understanding could also help in occupational health management, sue 


as reducing work-related musculo-skeletal disorders (MSDs) 


13, 


18 
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While the role of muscle fatigue in MSDs is not well understood, it has 
been established that when muscles fatigue the torque production strategy 
changes 19]. For instance, the body has to use muscle fibers that are more 
vulnerable to fatigue or use different muscle groups that are less efficient 
to maintain the level of torque demanded. If a fatiguing operation goes on 
without sufficient rest, such adaptations may increase the risk of MSDs. 

Muscle fatigue should, ideally, be predicted to aid in designing or evalu¬ 
ating repetitive industrial operations. Such a predictive model would allow 
optimizing working and resting phases, and could enhance long term health. 
Yet, such fatigue prediction models are difficult to develop. It is well known 
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that muscle fatigue is a complex phenomenon with multiple sources, of which 
the relative contributions remain unclear. For example, fatigue origins may 
be mechanical (peripheral fatigue: muscles are not able to respond to the 
input signal), neurological (central fatigue: the signal sent to the muscles is 
noisy or reduced), psychological, and/or other sources. As a consequence, 
predictive muscle fatigue models remain limited to restrictive assumptions 
and situations and focus on limited parts of the body. 

els have been previously proposed for muscle fatigue predic- 
. Most models, though, have focused on static situations: a 
desired torque at a fixed body configuration, which is an oversimplification of 
workplace tasks. In reality, the ability of torque production varies with body 
configuration due to the length-tension relationship. This paper proposes a 
more general muscle fatigue model, extending the static case to quasi-static 
situations, in which the subject’s posture can be changing, hence a vary¬ 
ing torque capability. It is worth noting that this new modeling approach 
cannot yet be considered as dynamic as the torque capability is assumed to 
depend primarily on the joint configuration alone. In dynamic situations, 
this capability depends on both joint configuration and joint motion (i.e. the 
force-velocity relationship in muscle force production). Such data will be 
required to establish a fully dynamic fatigue model in the future. 

To develop and analyze a quasi-static muscle fatigue model, we simulate 


arm motion during a repetitive flexion-ex 
task is common in industrial applications 


ension task. Even though this 


1] and daily life activities ll|, few 


studies have examined it with a varying body configuration 


16 


, 2, 3]. In this 


paper we first set the framework and assumptions for modeling. The arm 
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configuration data is then simulated and used to set up the model. Finally, 
the implication of the study findings is discussed. 


2 Problem setting for push/pull operations 


2.1 Definition of Push/Pull operations 

A push/pull operation generally defines both an action of pushing during 
which the force generated by the hand is oriented away from the body and 
an action of pulling during which the force is oriented toward the body. 
In daily life and in industrial environments, push/pull operations are very 
common, such as in objects displacement including up and down lifting or 
left and right shifting, drilling operations, buttons pushing, and so on. 


actions of pushing and/or pulling are an important source of MSD 


lese 


19] 


Indeed, NIOSH reported that about 20 % of overuse injuries were associated 
to push/pull operations [l?]]. Moreover, almost 8 % of back pain and 9 % o 


sprains and strains at the back are associated with push/pull operations 


18 


JLO]. In addition, heavy loads acting on the lumbar discs were detected when 


performing these activities 


19], which could also contribute to MSDs. The 


International Standard ISO 11228-2 


.6] recently proposed a first initiative to 


limit the professional required forces in push/pull operations with the aim of 
reducing or prevent MSD associated with this activity. 

Push/pull operations may be divided into three kinds, depending on the 
groups of muscles in which action is required to perform the motion: up- 
down; left-right or back-and-forth movements. In this paper, we will focus on 
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the back-and-forth movements, which corresponds to the majority of studied 
tasks in the literature and the most used push/pull operations in daily life. 

Another key component of a push/pull operation is the external load. 
While the whole body is involved in all push/pull operations, it is assumed 
here that, for the case of light external loads, the main groups of muscles are 
in the arm. Thus we focus on the arm movement with a light load push-pull 
task, while the movements of the rest of the body are neglected. 

2.2 Working assumptions 

We will set some working assumptions to simplify the study and address the 
problem of fatigue modeling in the frame of repetitive drilling operations. 
The study is made from the point of view of external forces. Specifically, the 
local effects of each working muscle will not be addressed, and instead the 
effects at the joint level (shoulder and elbow) will be considered for pushing 
or pulling. Assumptions made include: 

A1 Only light external loads are involved, and only the arm muscles are 
activated when pushing and pulling. Other muscles of the body will be 
neglected. 

A2 The arm movement will be limited to two joints: the shoulder and the 
elbow. The forearm, the wrist joint, and the hand will be considered as 
a single rigid body. This assumption is considered valid for negligible 
movements of the wrist, such as in drilling operations using back and 
forth movements. 
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A3 The group of muscles used for pushing will be considered inactive when 
the group of muscles used for pulling is active and vice versa. This 
assumption is equivalent to ignoring muscle co-contraction, or more 
precisely the use of antagonistic muscles does not contribute to their 
fatigue. 

The push-pull task used here is composed of two distinct phases. The 
pushing phase is generated by two groups of muscles: one for shoulder flex¬ 
ion (deltoids) and one for elbow extension (triceps). The pulling phase is 
generated by the two groups of muscles: one for shoulder extension (sub 
scapula) and one for elbow flexion (biceps). The model of fatigue will be 
applied to each active group of muscles. 

2.3 Cyclic motion 

A cyclic motion involves a repetitive movement where the beginning and 
ending positions, velocities and accelerations are identical. Let us consider 
a repetitive push/pull cycle with a regular periodicity noted T, such as the 
one illustrated in Fig. Q] Let T push and T pu11 denote respectively durations of 
the pushing and pulling phases durations with T = T push + T pu11 . 

Let i denote the number of completed cycles of the repetitive motion 
since the beginning of the work, and let t push and t pu11 denote the overall 
times counted from the beginning of the work respectively for pushing and 
for pulling. At the end of the work, we will have tT — tf — to where to and tf 
are respectively the starting and ending instants. If the work started with a 
pushing phase such as illustrated in Fig. Q3 then at the end of each pushing 
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Figure 1: Representation of the group muscles activity during push/pull 
operations 


phase we have f push = (£ + x^P ush an d at the end of each pulling phase we 


have f pu11 = fT pu11 . 

For an instant t E [t 0l tf], we have t = t push + t pu11 . Two situations may 
occur, depending if t is in a pushing or a pulling phase. If t is in a pushing 
phase, then f push = t — t pu11 = t — £T pul1 . The corresponding function for t push 
may be expressed using £ = (tf — t 0 )/T. 



( 1 ) 


If t is in a pulling phase, then we have f pu11 — t — t push = t — {£ + 1) • T push , 
which may be expressed with the following form. 



( 2 ) 
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3 Muscle fatigue model for push/pull opera¬ 


tions 


3.1 Formulation of general muscle fatigue behavior 


Using the assumptions in section 12.21 each group of muscles and its fatigue 
may be studied and its fatigue evaluated separately from the other groups. 


The fatigue model proposed in a previous study 
as the following differential equation. 


13] is adopted and defined 


dT cem _ Tj oint 


dt 


MVC 


•r r 


(3) 


where k denotes a constant fatigue parameter and T cem is the current maximal 
capacity for a group of muscles to generate a joint torque, and is considered 
as a characteristic value to evaluate the muscles fatigue. Tmvc (Maximal 
Voluntary Contraction ) is the initial value of T cem , measured at the beginning 
of the experiment. Its value depends on the operator’s anthropometry and 
body configuration as will be detailed in section 13.21 

Tjoint denotes the desired torque vector at the considered joint. It is 
composed of a portion coming from the body movements T* and a portion 
coming from the external load Tj )ext introduced as separate variables: 


Tjoint Tj T Tj^ext % 1 , . . . , 77 . (4) 

with n being the number of joints observed in the motion (n = 2 when only 
considering the s). T,; (i — I.. 77 ) are calculated from the values of 6, 6 and 6 






using the Lagrangian formalism [9] 


r, = 


d dL dL 
dt d6i d6i 


i = 1, 


,n 


(5) 


where L = E — U is the Lagrangian, E being the kinetic energy and U being 
the potential energy of the complete system (see 15] for the development ap¬ 


plied to the current problem). r ?iext is directly related to the desired external 
loads grouped in M load . Then r joint = r joint (0, 0, 0, M load ). 

Equation (j2J) may be re-written in the following form. 


r c em (t) = _ k ^ r joint (0(t) , 0(t), 0(t), Mload) 


i(*) 


MVC 


(9(f)) 


( 6 ) 


7 rt Rjoint(^( tjl )’^( t/l )’^( n )’-^load) ^ 

r cem (t) = C ■ e~ *o r Mvc( 9 W) u 


where C is an integration constant. 


(7) 


3.2 Formulation of joint capacities 

One of the key parameters of the fatigue model (0) is the maximal torque 
Tmvo I 11 static situations, in which a force is applied by the operator onto 
the environment with no body movement, the value of Fmvc varies according 
to its current body configuration. A predictive model to calculate Tmvc w as 
proposed by Chaffin et al. 4] and is summarized in Table [T] which gives the 
predictive functions for Tmvc at the elbow (index e) and shoulder (index s ) 
joints. Two parameters are considered as inputs: the joint angle value 9 and 
an adjustment gain G that accounts for the gender of the operator. Several 
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observations may be made based on the functions in Table [Q 


1. Joint capacity behaves differently for flexion and for extension move¬ 
ments; 

2. Predictive models for flexion are more complex than for extension; 

3. Only the shoulder extension model does not depend on the configura¬ 
tion of another joint. 


Joint / movement 

Joint capacity (r MVC ) 

G 

Male 

Female 

Elbow flexion 

(336.29 + 1.5446b - O.OO850 e 2 - O.50 S )G 

0.1913 

0.1005 

Elbow extension 

(264.153 + 0.5756b - O.4250 S )G 

0.2126 

0.1153 

Shoulder flexion 

(227.338 + 0.5256b - O.2960 fl )G 

0.2854 

0.1495 

Shoulder extension 

(204.562 + 0.0996b) G 

0.4957 

0.2485 


Table 1: Chaffin’s model for joint capacity in static situations [I] 


The model inputs are the measured shoulder and elbow constant joint 
angles, with the value 0 matching a straight arm in a classical standing 
configuration. From the joint angles, the respective shoulder and elbow re¬ 
spective capacities can be estimated, assuming this motion is relatively slow. 
The obtained joint capacity values will then be used for the fatigue model (JT]). 

It is worth noting that static strength data are being used here to model 
a dynamic task. While some dynamic strength data have been reported 
previously, such as in Frey et al. ?J, comprehensive dynamic strength data 
for the shoulder are not available. For simplicity, we used the same approach 
for both the elbow and the shoulder joints, based on Chaffin’s predictive 
functions. As a consequence, the proposed model is considered “quasi-static”. 
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3.3 Formulation of the push/pull fatigue in static sit¬ 
uations 

The model dTJ) is applied separately for each muscles group, for pushing and 
pulling. 


rt^ ppush^ 
k rt pu11 ppull , 

r?„m(i pu “) = r pu 5 c -e~ ? sf£ /o 


( 8 ) 


Because of Assumption A3, the model of fatigue of the arm for a push/pull 
operation is represented by piecewise continuous functions, in which each 
muscle group has increasing fatigue when working and remains at the same 
level of fatigue when relaxing (Fig. [2]). Recovery is not taken into account in 
this initial model development, which means that the subject should be less 
tired than predicted by the proposed model. Still, this model of fatigue leads 
to major observations on how human beings deal with fatigue independently 
from recovery. Depending on whether the current time t is in a pushing or a 
pulling phase, the model takes the two cases into account as follows. 


r c em (t) = 


r?em(^ push ), ^ £ pushing phase 
P cem (t pu11 ), t E pulling phase 


(9) 


When substituting (JT|) and (|2J) into (P), the following model is obtained. 


r c em (t) = 


ypusli / f 


rp'iH l + 

pom 1 ^ 


t S~ t o 

T 

t f~ t o 


T pu11 ) , t E pushing phase 


+ 1 


T push j , t E pulling phase 


( 10 ) 
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Figure 2: Schematic representation of fatigue in pushing and pulling groups 
of muscles (r cem ) at a joint during a repetitive push/pull operation in quasi¬ 
static situations 

It is important to note here that the two functions T^ 11 and are 
continuous. On the contrary, r cem (t) is not a continuous function as illus¬ 
trated in Fig. [21 which illustrates the model of fatigue (jTU[) and the previous 
comments for one joint animated by one pushing and one pulling groups of 
muscles during a repetitive push/pull task. Starting with the pushing phase, 
only the group of muscles dedicated to pushing is activated, while the pulling 
one is at rest (constant value of fatigue characterized by r cem ). When in the 
pulling phase, only the group of muscles dedicated to pulling is at work, and 
the pushing one is at rest. 

3.4 Formulation of push/pull fatigue in quasi-static 
situations 

Quasi-static situations denote the case when the arm is in motion and Fmvc 
varies with joint angles as defined in the section 3.2. These values do not 
remain constant but evolve with time together with the evolution of the 
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arm joint angles: r MVC = r MVC (0(t)). r MV c cannot be taken out of the 
integral in (J7|). In this case, the integration constant C is set by considering 
WO) = TMVC (0o) where 0 O = #(to)- The fatigue model in quasi-static 
situations is then expressed as follows. 


Tcem(^) — TmVc(0C 


, rt r joint( 0 ( n )>^( u )^( u )> M load) , 
g r MVC( 0 ( n )) U 


( 11 ) 


Similarly to the previous static case, the quasi-static fatigue model (TTT1) is 
applied successively to the pushing and pulling groups of muscles during the 
push/pull operation. Recovery during the resting phase is not considered. 


ppush^push^ 


r E“m( fP ““) 


CcW ■ e 

r£Sc(«») ■ <= 


fc/ 0 tPUSh 


ppush 

joint 


(S(u),e(u),S(u),M load ) 

r MVc( e (")) 


du 


fc/o PUU 


ppull 

joint 


(e(u),e(u),e(u),M load ) 


F puii 
1 MVC 


(»(«)) 


du 


( 12 ) 


Then the fatigue model for push/pull operations is obtained in quasi¬ 
static situations, similarly to static situations. 


i(f) = 



ppull 



t f~ t 0 

T 

trzL 0 

T 


ypull 


+ 1 


ypush^ 


t G pushing phase 
t G pulling phase 
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4 Application to push/pull tasks with back 


and forth arm movements 


4.1 Description of the push/pull task 

The tasks were chosen according to Ma’s work 14[, in order to compare the 
proposed quasi-static model to the static one. The arm of a male opera¬ 
tor of stature 188 [cm] and body mass 90 kg was modeled to simulate the 
performance of a repetitive push/pull task. The model of the task included 
the use of a tool positioned at the extremity of the arm and with mass of 
2 [kg]. The initial position of the hand was defined in the sagittal plane 
by P 0 = [P x ,o,Pz,o] [m] and the final position by P f = [P x j,P z j] [m]. The 
position {0, 0} corresponds to the fixed shoulder position (the origin of the 
reference frame). The operator generates a 20 [N] pushing force and a 10 [N] 
pulling force while tracking an horizontal line, back and forth. 

The push/pull cycle lasted 10 [s], with an equivalent time distribution for 
pushing and for pulling: T pus h = T pu n = 5 [s]. This kind of task represents 
an operation of classical horizontal drilling. It is illustrated in Fig. [3] for 
two different amplitudes of hand horizontal displacements: either with an 


amplitude of 20 [cm] (Fig. 3(a)), or an amplitude of 10 [cm] (Fig. 3(b)). In 
this figure, the left red circle at (0, 0) represents the fixed shoulder joint; the 
two segments, arm and forearm, are represented by the blue solid lines linked 
at the elbow joint represented by a moving red circle. The hand is at the 
end of the kinematic chain, represented by blue circles. The trajectories of 


the elbow are different for the two tasks. We can see in Fig. 3(a) that a 0.4 
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to 0.6 [m] horizontal displacement of the hand generates an elbow elevation 
of 10 [cm] and advancement of 7 [cm], while the 0.3 to 0.4 [m] horizontal 


displacement of the hand illustrated in Fig. 3(b) generates an elbow elevation 


of 3 [cm] and advancement of 3 [cm]. 



(a) The hand has a 20 cm amplitude 


(b) The hand has a 10 cm amplitude 


Figure 3: Arm movement during a drilling operation. T pus h = T pu n = 5 [s], 
= 20 [N], MC = 10 [N] 


4.2 Modeling of the operator’s arm 


The geometric model of the arm is summarized by the modified Denavit- 
Hartenberg parameters 5] given in Table [2j matching the kinematics repre¬ 
sentation in Fig. [H The motion only occurs in the sagittal plane (xo, zo), so 
the shoulder and elbow joints are represented by simple rotational joints and 
all the rotations angles are calculated along the yo lateral axis. 


Joint 

a 

d 

a 

r 

e 

Shoulder 

0 

0 

tt/2 

0 

0s 

Elbow 

0 

hu 

0 

0 

0e 


Table 2: DH parameters of the operator’s arm 
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Figure 4: Kinematic representation of the operator’s arm 

4.3 Trajectory generation 

The hand trajectory in the task space was described using polynomial inter¬ 
polation. The hand velocity and acceleration were both null at the beginning 
and end of each pushing and pulling phases. Beginning and ending positions 
denoted respectively by Pq and Pf. We obtain the following trajectory. 

P(t) = P 0 + p(t) ■ (P f -P 0 ), t 0 <t<tf 

with p(t) = 3 (t/tf) 2 — 2 (t/tf) 3 . The joint angles, angular velocities and 
angular accelerations were obtained from the hand trajectory profile using 
the inverse kinematics of the arm, leading to the joint profiles illustrated in 
Fig. [5] for one push/pull cycle. 

It is worth noting that defining the kinematics in the joint space is much 
easier when applying the model to a real human arm, as the center of rota¬ 
tion of the human arm may be difficult to locate and is needed to set the 
inverse kinematic model of the human arm. Nevertheless, the joints pro¬ 
files obtained from the definition of the kinematics in the task space offer 
more correspondence with the human real trajectories used to perform the 
described drilling operation, which means that defining the kinematics in the 
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(a) Elbow (b) Shoulder 

Figure 5: Shoulder and elbow joint angles, angular velocities and angular 
accelerations over time using trajectory generation in task space. The first 
phase is a pushing phase 


task space is more realistic. This kinematics description will be used in what 
follows. 

From the described kinematics, the joint torques at the shoulder, T s , 
and at the elbow, F e , can be calculated. Their respective evolution over 
time is illustrated in Fig. [6] for one push/pull cycle. The motion is set for 
an horizontal displacement of the operator’s hand from P 0 = [0.4, 0.1] to 


Pf = [0.6, 0.1] [m] which matches the case illustrated in Fig. 3(a] 



Figure 6: Shoulder and elbow joint torques rj oint for one operating cycle. 
The first phase is a pushing phase 
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Several observations may be made from Fig. [HI First, the torque at the 
shoulder joint is always greater than the one at the elbow. Second, the two 
curves show a discontinuity at 5 [s] due to the discontinuity in external loads 
between the pushing (20 [N]) and pulling (10 [N]) phases. During a given 
phase, the torques remain continuous but not linear: their values depend on 
the operator’s arm configuration. 


4.4 Evolution of joint capacity Tmvc 

As mentioned in section 13.21 Chaffin’s work 4] allows for determining the 
constant joint capacity r MV c for the shoulder and the elbow according to 
the arm configuration and depending on the type of desired external load 
(pushing or pulling). We have dissociated two working phases for a complete 
push/pull cycle: the pushing phase, in which the shoulder is in flexion and 
the elbow in extension, and the pulling phase: shoulder in extension, elbow 
in flexion. The evolution of the shoulder and elbow r M vc (Q a (t),0 e (t)) for a 
complete push/pull operation is illustrated in Fig. [0 



Figure 7: Evolution of Fmvc(^) at the shoulder and elbow joints for a 
push/pull cycle 
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Once again, a discontinuity is observed at the change between pushing 
and pulling phases. In this case, the discontinuity results because in our 
mathematical model only one T M vc(t) was needed to formulate push and pull 
(not to be confused with the individual muscle r MV c, which is continuous). 


4.5 Muscle fatigue 


The fatigue parameter k is assumed constant for a given operator and at 
a given joint, for any per 


from a previous study 
^elbow 0.24 


'ormed motion. The following mean values for k, 
14 were used in this framework: k s h ou ider = 0.17 and 
15j. This setting of k was obtained by an identification process, 
and was based on anthropomorphic data, maximal torques in flexion and in 
extension, and body dynamics representation. 

Using the expressions of the shoulder and elbow joint torques extracted 
from Eq. (0) and the matching joint capacities T M vc(^), the fatigue model 
can be applied to these two joints for a complete push/pull cycle. Figured] il¬ 
lustrates the resulting fatigue predicted at the elbow joint for three push/pull 
cycles. As the external load remains low, the decrease in joint torque capac¬ 
ity TMvc(t) is very slow and may not be easily observed in Fig. 0a), so a 
zoom-in on the pulling phases was realized and is shown in Fig. 0(b). Here, 
the increase in fatigue can be seen. Similar results were obtained for the 
pushing phases. 

From Fig. 0 we can see that the evolution of the fatigue (i.e., decrease 
in joint torque capacity r cem ) is not continuous at the transition between 
pushing and pulling phases. The main reason for this discontinuity is that a 
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Times (minutes) Times (minutes) 

(a) pushing and pulling phases (b) Zoom-in on pulling phases 

Figure 8: Simulation of the decrease of muscle capacity r cem at the elbow 
joint during repetitive push/pull operations for three cycles. The first phase 
is a pushing phase 

single r MV c was formulated for push and pull muscle groups. 

The evolution of the fatigue within a pushing or a pulling phase mainly 
depends on the values taken by the torque Tj oin1: . The resulting behavior is 
then different for a pushing phase or for a pulling phase as Tj 0 i nt is different 
according to the considered phase. 

The effect of fatigue for such light loads should be evident after a long 
period of push/pull operations. The model was thus simulated for 37 minutes 
of exercise, and Figure E3 illustrates the resulting fatigue predicted at the 
shoulder and elbow joints. For these two graphs, the bottom curve represents 
the oscillations of the desired torque Fj 0 i nt , which remains always the same, 
and the top curve represents the oscillations of the available torque r cem , 
which decreases with the increasing muscle fatigue. For both joints, the 
fatigue in pushing operation was more important than during the pulling 
operation. This is explained both by a less important joint capacity r cem 
and a higher external load for the pushing phases. Using such graphs to 
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predict potential MSD risks, such risk would seem most important when the 
two parts join (i.e. at around 26 minutes for the elbow joint and at around 
35 minutes for the shoulder joint in this case). In this example too, MSD risk 
would seem higher at the elbow joint. Still, as the recovery model is not yet 
included in this approach, and thus in reality any risks would likely involve 
longer periods of activity. 



0 5 10 15 20 25 30 35 40 0 5 10 15 20 25 30 35 40 

Times (minutes) Times (minutes) 

(a) Elbow (b) Shoulder 

Figure 9: Simulation of the decrease of muscle capacity at the shoulder and 
elbow joints during repetitive push/pull operations, and with the hand mov¬ 
ing from P 0 = [0.4,0.1] to Pf = [0.6,0.1]. The bottom curve represents 
the oscillations of the desired torque rj oint and the top curve represents the 
oscillations of the available torque r cem with respect to the push and pull 
actions 


4.6 Discussion 

4.6.1 Effect of constant or variable joint capacity 

In this paper, the joint capacity r M vc was introduced as a variable to build 
a model of fatigue in repetitive push/pull operations. To highlight the differ¬ 
ence between the use of T M vc as constant or variable (i.e. depending on the 
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posture), we compare the fatigue graphs obtained in the two situations. The 
minimum value of Tmvc obtained for all the simulated configurations was set 
as the reference constant value; this value will produce the greatest muscle 
fatigue. Figure ITOl shows the level of fatigue in the situation when r MVC is 
set as a constant. 



(a) Elbow 


(b) Shoulder 


Figure 10: Simulation of the decrease in muscular capacity of the elbow and 
shoulder when Fmvc is a constant. The bottom curve represents the oscilla¬ 
tions of the desired torque Fj oint and the top curve represents the oscillations 
of the available torque r cem with respect to the push and pull actions 


The difference between the curves obtained by considering a variable value 
(figure E} and a constant value (figure [TUI) is very small in this framework. The 
level of fatigue is slightly greater in the case when Fmvc is constant, because 
the worst situation was considered. An increased risk at the elbow appears 
after 23 minutes instead of 25 minutes, which is almost identical. So, it 
seems that the first approximation of a constant value for the joint capacity 
is adequate to evaluate the muscle fatigue, if including a slight allowance 
for the prediction of MSD. However, the current work setting may have an 
influence on this observation, and a general sensitivity analysis should be 
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conducted to establish a more general conclusion that using a constant T M vc 
is adequate to evaluate the muscle fatigue. 


4.6.2 Effect of two push/pull tasks 

The hand movement amplitude was previously considered to be 20 [cm], 
moving from Pq = [0.4, 0.1] to Pf = [0.6, 0.1] [m]. If the initial and final 
positions are modified to P 0 = [0.3, 0.1] and Pf = [0.4, 0.1] [m], as depicted 


in Fig. 3(b), the resulting fatigue curves are shown in Fig. |TT] 



0 10 20 30 40 50 60 0 10 20 30 40 50 60 

Times (minutes) Times (minutes) 


(a) Elbow (b) Shoulder 

Figure 11: Evolution of fatigue for the hand moving from P 0 = [0.3, 0.1] to 
Pf = [0.4, 0.1]. The bottom curve represents the oscillations of the desired 
torque rj oint and the top curve is the oscillations of the available torque r cem 
with respect to the push and pull actions 


As expected, the two push/pull task configurations generate different fa¬ 
tigue. The simulation was run this time for a sixty minutes of exercise. Still, 


Fig. 11(a) demonstrates that the fatigue curve for the elbow in the pulling 
phase is slightly increased: this may be observed on the upper values of the 


top oscillations (Fig. 11(a)), starting at 71 [m.s] and ending at 73 [m.s]. This 
may seem an unrealistic result, because performing an operation should in 
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theory reduce the joint capacity r cem . To better understand the reason for 
this behavior, let us observe the torque evolution rj oint and joint capacity 
Tmvc f° r the studied operation, illustrated in Fig. [T2J The second phase of 
Figurdl2l is the pushing phase. Figurdl2fbi demonstrates that the torque 
level is close to zero at the elbow joint, thus minimal fatigue is generated for 
the push muscle group at the elbow joint. This particular finding may have 


- shoulder MVC 
elbow MVC 


4 6 

Times (seconds) 

(a) r M vc 



(b) Feint 


Figure 12: Evolution of the input torques at the elbow and shoulder for the 
hand moving from P 0 = [0.3, 0.1] to Pf = [0.4, 0.1] 


workplace implications. While the differences between the two push/pull 
task configurations used in the simulation are seemingly small, the end effect 
on fatigue can be dramatic (i.e. for the push muscle group at the elbow 
joint). Thus by configuring workplace tasks carefully, it may be possible 
to reduce the workload and/or fatigue for the most vulnerable muscles and 
reduce MSDs. 
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5 Conclusions 


In this article, a quasi-static model of muscle fatigue was proposed as an 
extension of a previous model of static fatigue, by incorporating a variable 
joint capacity r MV c as a function of the operator’s posture. The model was 
applied to simulate repetitive push/pull operations with light external loads 
that may be observed in an industrial framework, such as a classical drilling 
operation. From the simulation of fatigue during one task configuration, 
we found that fatigue of the elbow appears faster than the shoulder. This 
result is in agreement with the observation that most arm MSDs appear at 
the elbow joint. Indeed, the muscle strength of the elbow is lower than at 
the shoulder. However, for the second task configuration, we found minimal 
fatigue for the elbow. This last observation demonstrates that it is possible 
to reduce MSDs by optimizing workplace tasks. 
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